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We offer a general treatment of electromagnetically induced transparency (EIT) in a five level 
system consisting of four metastable ground states. Two additional RF/microwave fields coherently 
couple the two ground states of the standard EIT A atom to a pair of additional hyperfine states 
in the ground state manifold, generating two sets of dressed states that interact via the probe and 
control lasers, which couple to the electronic excited state. These new hyperfine dressed states 
manifest themselves in the linear optical susceptibility of the probe as new resonances in addition 
to the Autler-Townes doublet characteristic of EIT. In particular, we show that the existence of 
two new narrow resonances, whose width are limited only by ground state decoherence, appear 
inside the normal EIT transparency window. We show that by controlling the intensity of these 
RF/microwave fields, one can engineer both the position and width of these narrow resonances and 
thereby exercise additional control over both the dispersion and group velocity of the probe. 



I. INTRODUCTION 

Electromagnetically induced transparency (EIT) [l[ is 
a quantum interference effect in which coherence be- 
tween atomic states is used to reduce the absorption in 
a window around an atomic resonance, while simultane- 
ously generating large dispersion and third order nonlin- 
ear susceptibilities within the induced transparency win- 
dow [3, Q. The simplest realization of an EIT system 
consists of a driven A-type atomic system wherein the 
excited state is coupled to one ground state — the auxil- 
iary state — via a a strong control beam; meanwhile the 
system is probed by a weak laser beam near resonance 
with the transition between the excited state and the as- 
yet uncoupled ground state. 

In these simple systems, the characteristic EIT inter- 
ference effects arise because the control beam induces a 
splitting of the excited state into a doublet of "dressed" 
states, providing two interfering excitation pathways for 
the probe beam. The null in the probe absorption at the 
normal atomic resonance frequency is a result of destruc- 
tive quantum interference between these two pathways. 
The resulting dispersion has very large slope for frequen- 
cies within the transparency window around the absorp- 
tion null, leading to extremely slow probe beam group 
velocities [4, 5]. Slow light propagation of this sort has 
been observed in a variety of media, including hot atomic 
gases p and atomic Bose-Einstein condensates [7|- 

Several groups have expanded on the basic EIT system 
to further explore the optical properties of coherently pre- 
pared atoms. In general, these modifications to standard 
EIT have involved coupling additional states to the three 
levels of standard EIT, via new lasers or RF/microwave 
fields. One class of these have led to EIT's sister phe- 
nomenon, electromagnetically induced absorption, which 



occurs when the primary ground state is strongly cou- 
pled to an additional, off-manifold state, forming an N- 
configuration [3, Q. Others involve levels on the same 
electronic manifold as the EIT ground states. One such 
modification considers a four level configuration, in which 
an additional ground state is coupled to the auxiliary 
state via an RF or optical field. The resulting doublet 
constructively interferes with the excited state to produce 
a new absorption feature in the EIT spectrum, with prop- 
erties controllable by tuning the Rabi frequency of the ad- 
ditional field [IJJ] • These so-called "interacting dark reso- 
nances" have been observed experimentally in several ex- 
periments [ll], Q2, 1 and they have been suggested for 
the coherent control of group velocity in the UV regime 
[14j | . Additional theory and experimental work has shown 
that the feature arising from the interaction between the 
dark states persists, even in the presence of considerable 
Doppler broadening [T^, [l6| • 

More recently, in an attempt to explain experimental 
results [13| , the absorption spectrum of a similar four 
level atom was considered, where now the additional 
state was coupled to the ground state that is probed 
by the probe laser [17| . Other modifications have in- 
volved coupling the two grounds states of the A atom (or 
other 3-level atomic configuratio ns, such as cascades) to 
each other, via an RF field [13, Il3, l23|- ^ n these cases, 
a variety of new phenomena are possible, including effi- 
cient frequency conversion for the probe beam. Alterna- 
tively, rather than introduce additional hyperfine levels, 
recent work has explored the modifications to EIT result- 
ing from coherent tunneling of A atoms in a double-well 
Bose-Einstein condensate, forming an effective six-level 
configuration [2l|. Here, space- like separated atoms can 
be used to modify the optical response of a sample in one 
well. 



In the current contribution, we consider the more gen- 
eral coherent interactions that occur when both of the A 
atom ground states are coherently coupled to additional 
states. These interacting dressed states generate new fea- 
tures in the EIT absorption spectrum, whose widths and 
locations in frequency space can be controlled by inde- 
pendently modulating the Rabi frequencies of the two 
new coupling fields, providing markedly better control 
over the optical properties of an atomic sample than pos- 
sible with previous implementations. Moreover, we show 
how, by controlling the intensity of one of the new cou- 
pling fields, one can generate group velocities almost 100 
times slower than otherwise possible in a given EIT sys- 
tem with a fixed control laser intensity. Whereas other 
studies have considered replacing either one or the other 
ground state with a driven doublet, no general treatment 
has appeared in which both ground states are driven. 
Thus, the phenomenon of the interaction between the 
doublets — which is what permits the fine tuning of the 
resonance locations within the transparency window — 
has yet to be explored. The ability to tune the relative 
position of the narrow resonances is what allows for the 
enhanced control over the dispersion within the trans- 
parency window. The current paper includes an ana- 
lytic statement of the linear susceptibility of the five level 
atom, from which both the absorption coefficient and dis- 
persion can be simply derived. 

The remainder of this paper is organized as follows. 
In section UH we describe our model for the five level 
atom, dressed by two RF/microwave fields and a strong 
control laser. We derive the master equation in a par- 
tially dressed basis and show that in the steady state, 
the coherence term in the master equation relevant to 
the optical properties of the probe beam permits an an- 
alytic solution. In section IIII1 we derive the system's 
linear susceptibility, x , from which the optical proper- 
ties of the sample can be extracted. Here we also analyze 
the dispersion and group velocity of the probe laser. Fi- 
nally in section IV, we comment on present experimental 
constraints. An appendix includes the full equations of 
motion for the relevant terms of the master equation. 



II. THE MODEL 

We consider a five level atom as depicted in Fig. 
1 where two pairs of ground states {{\b}, \b'}} and 
{|c),|c')}) interact with a single excited state, \a). Di- 
rect transitions between these ground states are assumed 
to be electric dipole forbidden (ie. \b) </-> \c') or \b) ^A |c), 
etc.); moreover, to avoid possible degeneracies, we as- 
sume that any degeneracy of the states is lifted by an 
external magnetic field. In analogy to a standard EIT 
configuration, a strong laser of field amplitude £^ and 
of fixed frequency w M propagates near the \c) «-> \a) 
transition. Here, we study the propagation of a weak 
probe laser (£ p <C £ ^) of variable frequency io v near the 
\b) <-> | a) transition. The couplings between the atomic 




FIG. 1: Our 5-level model. \b), \b'), |c), and |c'} are all in 
the same ground state manifold; \a) is an excited state. In 
analogy to standard EIT, \a) and \c) are coupled by a strong 
control beam, Q^. \b) and \b') and \c) and |c') are coupled by 
RF/microwave fields, Qb and il c , respectively. A weak probe 
beam, fi p , propagates near the \b) <-> \a) transition. 



levels are moderated by their complex Rabi frequencies, 
which are defined in terms of the lasers driving the tran- 
sitions. The \c) <-> | a) Rabi frequency, 51 M , is given by 
M7 M e~ l< ^ = £fiD ac while \b) <-> |a) Rabi frequency is 
similarly given by h£t v e~ %l ^ v = £ p D a b. D$j = e(i\x-e\j) is 
the dipole moment matrix element in the direction of the 
laser polarization, e, for the \j) <-> \i) transition. Addi- 
tionally, two RF/microwave fields drive magnetic dipole 
transitions between members of each hyperfine ground 
state pair: one, of frequency u>b and Rabi frequency fib 
drives the |6) <-> \b') transition; another, of frequency u) c 
and Rabi frequency Q c drives the \c) «-> \c'). From here 
on we will refer to the new fields as RF fields for the sake 
of simplicity. f2 M , f2&, O c , and f2 p are all taken to be real. 

We take the base vectors for the five atomic levels \a), 
\b), \b'), |c), and \d) to form a basis for the relevant 
subspace of our overall Hubert space. Then, in a frame 
rotating at the frequencies of the fields, defined to avoid 
explicit time-dependence, we can write the state vector 
as 



|*(i)> = MtM+Mt)\b)+i>b>(t)\b')+Mt)\c)+ipc>(t)\c'), 

(1) 

where 



& = e^e - *"**^ 






-i<l> p +i<j> b —iv v t+iv h t 



h 



-i<t>u 



Vc 



ij) c , =e -^+^ e -^ t +^ t ^ c ,. 



(2a) 
(2b) 
(2c) 
(2d) 



The Hamiltonian in the rotating frame is then given by: in the {\b), \b')} plane, 9b'- 



U =- (u) a \a)(a\ + (w b + v P )\b)(b\ + {ojy +u p - u b )\b'){b'\ 
+ (u> c + ^)|c)(c| + (uic> + Uft - u c )\c'){c'\ - M |a)(c| 



n 6 |&')<6|-nc|c')(c|-np|o)(6|) + h.. 



(3) 



To accommodate the possibility of strong fields driving 
the hyperfine transitions, we will want to keep fib and fl c 
to all orders in our calculation. To do so, it is convenient 
to work in a partially dressed basis, in which the sub- 
spaces of the effective Hamiltonian corresponding to the 
two ground state doublets are diagonalized. We will ex- 
plicitly diagonalize the {\b), \b')} subspace; the {|c), |c')} 
subspace behaves identically mathematically and there- 
fore all results can be inferred from those for the {\b), \b')} 
subspace. The {\b), \b')} subspace of the Hamiltonian can 
be written as a sum of its diagonal and traceless parts: 



1~Lw — — ( {2uj b 



A h 



2i/ p )I - 



-A 6 -n b 



(4) 



where Aj, = u>b' — <+>b — ^b is the detuning of the RF field 
driving the \b) «-> \b') transition. 

The matrix diagonalizing Hbb' will be a member of 
SO(2), so it can be written in terms of a rotation angle 



D h = 



where 



cos 8b — 



cos Ub sin Ub 
- sin 9b cos Ob 



l + A b /flf^ :! 



(5) 



(6a) 



n^^f 


(6b) 


nf^JAl + nl 


(6c) 



Applying this transformation produces dressed states 
\B) = cos9 b \b) + sm6 b \b') and \B') = -sin6 b \b) + 
cos 9b | b') in the diagonalized basis with the energy eigen- 
values ±hfl e b ff /2. 

Combining the transformations for the two subspaces, 
Db and D c , into a block diagonal matrix, we find the di- 
agonalization matrix, D, which can be used to determine 
the partially dress Hamiltonian: 



J 



DHrf = - 
2 



( 2oj a 

— flp cos 9 b 
n p sin Ob 

-fl^cos9 c 
\ n^sinOc 



2io h 



-fl p cos 9b 



2v„ 



nf 



2u) h 



flp sin Ob 


A b + 2v p 





r 



Qf 



2lu c 



-flftCOSt 







2v,, - n 



oil 



2U)r 



Qft sin C 






A n + 2u u 



Of J 

(7) 



At this point we move from a wave function descrip- 
tion of the atom to a density matrix, py, that al- 
lows us to incorporate decay and decoherence into our 
model. Since the linear response of the probe is deter- 
mined by the coherence between the states |o) and \b), 
Pab = cos ObPaB — sm9bp~ a B', the primary task will be 
to calculate the steady state solutions of the partially 
dressed density matrix elements, p a B and p a B' to linear 
order in the probe field. 

In this paper we consider a closed system in which pop- 
ulation from the excited state \a) decays to the included 
ground states via spontaneous emission and the ground 
states are stable with respect to decay. Consequently, 
population is conserved within our 5 level basis and there- 
fore we do not include pumping of the states. Note this 
restricts our discussion to cold atomic gases where atoms 
do not quickly pass through the region of interaction with 
the lasers, which allows us to ignore pumping and de- 
cay due to atom transits through the interaction region. 
For the ground states, we include pure dephasing deco- 



herence between states due to collisions, stray magnetic 
fields, etc. 

The total spontaneous emission rate is given by 7 a and 
we assume that \a) can decay with equal likelihood to 
\b), \b'), or \c) at the rate 7 Q /3, but because of opti- 
cal dipole selection rules, it cannot decay to a fourth 
state. We denote pure dephasing decoherence rates be- 
tween two states by 7^ . We make the following assump- 
tion that the difference in the dephasing between \a) and 
both members of the doublet {|6), |6')} are much smaller 
than 7 a , i.e. \^ a b — 7o6'| *C 7a/6 so that we can set 
lab = lab' = la /6 + lab for the total decoherence rate. 
This approximation is justified by experimental param- 



eters where i a 



10' 



while collisional decoherences 



rates can be anywhere from a few tens of KHz down to 
hundreds of Hz or less for quantum degenerate atomic 
gases or with buffer gases [22|, [2J, |2J| . Additionally, we 
take the decoherence between the dressed state doublets 
to be of two types, which we allow to vary independently: 
we have i cb = lev = lc and j &b = Idv = lc>- 



These two assumptions lead to a considerable simplifi- 
cation of the equations of motion in the partially dressed 
basis with the relevant contributions being: 



PaB ~ —"labpaB 
PaB' ~ —labPaB' 

- (7c cos 2 1 



PCB 



- JC sin 2 9 c )pCB 



- (lc - lc) cos 9 C sin 9 c pcB 
Pcb ~ ~(lc sin 2 9 C + 7c cos 2 9 c )pc'B 

- (lc - lc) cos 9 C sin 9 c p CB 
Pcb 1 ~ -(7c cos 2 9 C + 7c sin 2 9 c )p C B' 

- (lc - lc) cos C sin 9 c pc>B> 
PCB' ~ -(7C sin 2 6> c + 7c' cos 2 9 c )p C 'B> 

- (lc - lc) cos 6> c sin 9 c pcB' 



thus, Eqs. [23] simplify to: 



,<9p 



Bfl 



«76b' 



dt 

.dpB'B' 



(pbb - Ps'B') 



llbb' 



dt 



(pbb - Ps'S') 



■ dpBB' , n ilw,- 

»-^-=(-n 5 --2-) PBB , 



Hi*' - 
—^-PB'B 



dt 



(C , *lbb> V - 



llbb' - 
—Z-PBB' 



Since dpBBJdt and dps'B' jdt depend only on the differ- 
ence between pss and pb'B' , these equations are solved 
by equal constants — since we have assumed that initially 
all of the atoms are in the {\b), \b')} manifold, we have 
Pbb = PB'B' = 1/2. As for the coherence terms, we can 
write their coupled equation of motion as a matrix: 



While this second assumption may seem strange, it per- 
mits both an analytic solution, and maintains sufficient 
nuance to allow exploration of the effects of different de- 
coherence pathways on the new resonances as will be 
shown later. Note that we do not include lac, lac' in 
this discussion since they do not enter into the solution 
for pab to first order in the pump field. 



We assume that all of the atoms are initially in ei- 
ther \b) or \b'). Since we are solving to linear order in 
the probe under the assumption Q p <C lab, many of the 
terms in the density matrix can be discarded since they 
only develop population at order (fl p /i a b) 2 or higher in 
perturbation theory. For example, p aa ~ ^p/lab ~ 
while {|c), |c')} manifold only develops population at or- 
der fi 2 f2 2 /7^ b which is also negligible. We can then take 
Pec ~ Pc'd ~ Pec' ~ Pac ~ Pad ~ 0, and likewise their 
complex conjugates. Note, however, that the assump- 
tion of a weak probe and thus of negligible occupation 
of \a), |c), and \d) amounts only to the assumption that 
il p <C lab- In particular, it places no constraints on the 
strength of the U and £l c . With these assumptions, we 
find two systems of equations that are decoupled from 
each other, which we include in the appendix. 



Also included in the equations in the appendix are 
those for the subspace \B), \B'), which we note are de- 
coupled from the other equations under the assumption 
that all terms higher than first order in f2 p are negligi- 
ble, implying that they can be solved separately from 
the other equations. The terms pbb, Pb'B', Pbb', and 
PB'B, which appear in the other equations, act as source 
terms. In the special case of an on-resonance field driving 
the | b) «-> \b') transition, we have 9 b — ir/4 and fi£ ff = fib', 



; d_ (PBB' 

'dt \PB'B 



Hbb' 



-n b 

i ibb' 

2 



'Hbb' 
2 

n b 



fpBB' 
\PB'B 



(8) 



The solution for pbb' and pb'B can be written as linear 
combinations of the eigenvectors of the matrix, which 
have the time dependen ce X\ = C\e { - lA ^ B ^ t and X2 — 
Q 2e {-iA-B)t with A = ^4(^)2.^/2 and B = lw / 2 . 

For t ^S> 766', both of these go to zero, which gives the 
steady state solutions pbb' — Pb'B = 0. 

To solve for p a B and p a B' , we assume that the 
{\B), \B')} manifold has been given time to reach the 
steady state we have described above. Plugging these 
into the master equation for p a B and p a B' , and assuming 
further that the RF coupling |c) <-> \c') transition is also 
on resonance, A c = 0, (giving 9 C — 7r/4), we find the 
simplified systems: 



.dpai 
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£l b . \ fL,V2 
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+ 4 PaB' 2 (7C " 
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(10c) 


where A p 


= L0 a - uj b - v p and A^ = 


- LO a - W c - ^ 


are the 



probe and control beam detunings, respectively. 

Eqs. [9] and [10] have identical structures; both can 
be solved analytically by writing the systems as matrix 
equations of the form dX/dt = —M.-X(t) + A and noting 
that equations of this form have steady state solutions 
limt^oo X(t) = M _1 • A. We then extract the steady 



state solutions to p a B and p a B' from the resulting vec- 
tors. 



III. OPTICAL PROPERTIES OF THE 5-LEVEL 

SYSTEM 

The complex linear susceptibility, expanded about the 
| a) ^> \b) transition, is given by 



(1) 2a(v)ND ab _ 
X = ? Pab- 



£o£p 



(11) 



X^ determines both the absorption coefficient, a(A p ) = 
kp^slx^iAp)], and the index of refraction, n(A p ) rs 
(1 + Kfx^^Ap)]) 1 ' 2 . Considerations arising from the 
particular experimental set-up would determine den- 
sity profile, <r(r), and number density of atoms, N, 
which are included to formally account for the partic- 
ular optical thickness of the sample. In our analy- 
sis we will focus on the reduced susceptibility x = 
eoftlabX^ /1D 2 ab Na(r) = j a bPab/^ P - We then arrive at 
the following analytic expression for x , 
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y(l) = 2^ 

X 2 



{Ap -A p + irtc ~ ^fc/2)(A M -A p + i lc - «&/2) - n 2 c /A 



+ 



(A p - i lab + fi 6 /2)((A„ ~A P + i~f C > - n t /2)(A M - A p + i lc - fib/2) - njj/4) + (A p - A p + i<y c > - fi b /2)ft2/4 
(A p - A p + v)c> + ^b/2)(A p - A p + i lc + Sl b /2) - n 2 c /A 



(Ap - i lab - fi 6 /2) ((A M - A p + i lc , + fi 6 /2)(A M - A p + i lc + Q b /2) - Q 2 /A) + (A A 



«7c 



n b /2)^l/A 
(12) 



r 



It is worth pointing out that in the limit that Q b , Q c — > 0, 
we recover the standard expression for the coherence for 
EIT in a A system, 



x (1) 



7a6 



llcb 



2 (Ap-ij ab )(A p -i% b )-(n il /2) 2 



Further analysis requires us to estimate the values for 
the important variables in the problem. The Rabi fre- 
quencies of the coupling laser and RF fields are experi- 
mentally tunable over a large range. For the spontaneous 
emission rate, we take 7 a = 10 7 s~ 1 and will measure the 
Rabi frequencies and detunings in units of j a b- For the 
ground state dephasing rates the range of values are lim- 
ited primarily by collision rates and therefore are tem- 
perature and density dependent. However for concrete- 
ness, we assume that 7 O0 , 7c > an d 7c are in the range 
10 3 - 10 4 s _1 . 

Examples of the real and imaginary parts of the sus- 
ceptibility are shown in Fig. [2] for A M = while Fig. 
[3] display the spectrum's dependence on the Rabi fre- 
quency fif,- We see that the presence of the additional 
levels manifest themselves as two narrow resonances lo- 



cated inside of the EIT transparency window. In general 
for arbitrary A by the new resonances are symmetrically 
located about A p = at the locations A p = ±fi£ ff /2. For 
fi^t,7a6 S> Q(,,fi c , 7c,7c, their shape is approximately 
Lorentzian, given by (for A M = 0): 



3[X 



o(i)l 



7afe^c 

2n u 



n 2 jn 2 + lc >hab 



(A p T n b /2) 2 + ( lab (n 2 c /nl + IC'/labW 

in the vicinity A p « ±Q b /2. In the limit as 7c — > 0, 
this expression reduces to a Lorentzian with FWHM of 
2 r y ab il 2 /n 2 , which agrees with what has been obtained by 
Lukin et al. [10( and Mahmoudi et al. [141 ] for interacting 
dark resonances in a 4 level system with no ground state 
dephasing. In the case of nonzero 7c, the widths of the 
features, 



r, 



2 lab (Q 2 c /nl + jchab), 



(14) 



is the sum of the 'power broadening' term 2^f ab tt 2 /£l 2 
and the dephasing rate for \c') while the height is given 
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FIG. 2: (Color Online) (a) The full spectrum of the suscepti- 
bility showing the imaginary part, S£[x ]> (blue dashed line) 
and real part 5R[x ], (red solid line), (b) A close up on one 
of the new narrow features. Here, fib = fi c = "fab/10 and 
Q ab — 2j ab , while jc = 7c = 0. 



by 



»[x (1) (±^ 
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2(7„6fi»+I#7C')' 



(15) 



These new resonances have a simple interpretation 
in terms of the dressed states \B) and \B'), which 
are coupled via the probe laser to the 3 level ststem 
\d) «-» \c) «-+ \a). The energies of \B) and \B') arc 
fux>B,B' = Swb ± ^b /2 so that in the absence of the con- 
trol laser, the \b) <-> \a) absorption line would be split 
into an Autler-Townes doublet located at uj a — ujb and 
Lu a — lob 1 - The \d) «-> \c) <-> |a) is isomorphic to a A atom. 
Again assuming A c = and A p = 0, we then have the 
following eigenstates for the {\a), \c), \c')} subsystem, 



\a+) = -=(smO\a) + \c)+ cos 9\c')) (16) 

v2 

|o_) = -p(sin0|a)-|c) + cos0|c / )) (17) 

|o ) = cos(9|a) - sin6»|c') (18) 

where tan 9 = —fl^/il c . The energies of the states \a±) 
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are _E± = ?iaj a ± hJfl 2 + 0§/2 while |ao) has energy 
£b = huj a . As one can see, |ao) is the same type of 
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FIG. 3: Imaginary part of the susceptibility, Q[x ] plotted as 
a function of A p and fi;,. The separation of the new features 
varies with fi;,, the Rabi frequency of the RF field coupling 
the states \b) and \b'). Here, we show the dependence for 
< fi b < j a b/10, with 7c = 7c = lab x 10" 3 , fi M = 2j ab , 
and fi c — jab/10- In the intensity plot, lighter color denotes 
larger Q[x ]■ 



dark state that appears in STIRAP and coherent popu- 
lation trapping. In this case, this induced dark state is 
a superposition of \a) and \d) but not \c) . Since it is 
decoupled from the control laser, there will not be any 
destructive quantum interference in the probe absorption 
for transitions to |oo). Transitions from the {\B), \B'}} 
manifold to |ao) will then exhibit absorption resonances 
at uj a — Ub,b'i which correspond to the new narrow res- 
onances. By contrast, destructive interference created 
by the control would lead to nulls in the absorption at 
Lu a — ojb,b' when Q c = 0. Fig. [4] shows a schematic 
diagram of the energy levels of the dressed ground state 
manifold {\B), \B')}, which are coupled to all three states 
of the excited state manifold {|a+), |a_), |oo)} via the 
probe. All in all there are six transitions that should 
appear as resonances in the absorption spectrum. The 
transitions to the two 'bright' states \a±) correspond to 
the main absorption peaks located at A p s=s ±f2^/2 for 
il M ^> O c , £!{,. Notice that each of these resonances actu- 
ally consist of a pair of resonances separated by a distance 
fib but only when fib > j a b can these pairs be individu- 
ally resolved as shown in Fig. [5] 

The independence of \ao) from |c) explains why the line 
widths of the new features only depend on the dephasing 
between states in the {\b), \b')} manifold and \a) and \d) 
but is independent of dephasing relative to state \c) as 
one can see from Eq. 1131 which is independent of jc ■ Fig- 
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FIG. 4: Energy level diagram that indicates transitions in- 
duced by the probe laser between the ground state manifold 
{\B), \B'}} and the excited state manifold {|a+), \a,~), |ao)}- 
Transitions to the dark state \ao) are indicated by dashed 
lines. The energy of the bare state \b) is also shown for refer- 
ence. 
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FIG. 6: (Color Online) Imaginary part of the susceptibility, 
Sj[x ], as a function of A p showing a close up of the narrow 
resonances. Here we show the dependence of the new features 
on the dephasing, <fc'- The broad dashes represent the case 
where jc' — 0, the short dashed line takes jq/ — 7 a (, x 10~ 3 , 
and the dotted line represents jc' = lab X 10 -2 . In all cases, 
tt b = Q c = 7ab/10 and Q M = 27 ab . 
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FIG. 5: Imaginary part of the susceptibility, Ss[x ] 
function of A p . Here, fi M = 2"/ a b, Qb = 2.27 a t, Q c = 
have been chosen so that all six absorption resonances are 
simultaneously visible. See text and Fig. [4] 
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[5] shows how S[%^] varies with jc 1 - Moreover, that the 
height of the features decreases with increasing 7c , as 
in Eq. Q1)J reflects the fact that the dephasing decreases 
the probability of transitions to the dark state |ao). 

These new narrow resonances offer the possibility of 
additional control of the dispersion, d$t[x^]/duj p , and 
group velocity, 

Vg{U p 



n+(uj p /2n)(dn[x^}/du; p ) 



inside of the normal EIT transparency window because of 
the large normal dispersion in the vicinity of the narrow 
resonances, whose positions and widths are controlled by 
£lb and f2 c , respectively. By controlling the RF Rabi fre- 



quencies, the narrow resonances can be selectively posi- 
tioned inside of the normal transparency window in order 
to locally control the dispersion over a relatively small 
frequency range. 

To be more specific, we note that in 'standard EIT' 
(note that 'standard EIT' refers to the case where 
Qf, = fi c = but with all other parameters, includ- 
ing the control laser, being the same), the observable 
effects of EIT including a transparency window and slow 
light occur when |fi M | 2 3> lablcb but at the same time 
d$l{x {1) }/duj p oc 2 7afc /|f7 M | 2 for A p « [2]. This implies 
that although lowering f2 M will decrease the group veloc- 
ity, the decreasing of the intensity of the control laser to 
l^^l 2 < labia, will result in a complete loss of the trans- 
parency window. Let us return now to our 5 level system 
and focus specifically on the case where f2 M , 7 a {, ^ ^h, T„ 
so that the narrow resonances are clearly visible in be- 
tween the Autler-Townes doublet induced by the control 
laser. In this case the narrow resonances have line widths 
given by Eq. 1141 which requires that Qb 3> T n in order 
to have a clearly defined window around A p = 0. The 
natural line width of the excited state only contributes 
through the power broadening term 2(fl c /n tJi )' 2 j a b <C j a b 
for Q c <C flfi and can easily be made much smaller than 
the ground state dephasing rates. As a result the primary 
limit on the dispersion and absorption in the middle of 
the resonances arises from the ground state dephasing 
alone, requiring Qb 3> ic ■ This limit is a much weaker 
condition than for the control laser in standard EIT when 
lab 3> leb and allows for much narrower transparency 
windows embedded inside of the standard EIT window. 

In our system, at A p = and assuming 
fib, ri c ,7c, lc ^ ^nilabi the absorption and dispersion 
are, to lowest nonvanishing orders of Qb, ^c, lc , and lc -, 
given by, 



is 



[x (1) (0)] 



lab 



&M 



and 



4 7afc (» 2 -»g) 2 + ( 7c ^+ 7c ^ 2 )» 2 

- nir - mi + (7ab7c , + ni)nmi + 



W 



a3?[x (1) ]/5A p - 7ah 



-87ab(7C^fo + {lc 



•37C')W 



7c fi; 



4 ^ 2 



ng(^ + n>)nj 



fi2 (8 7o6 ( 7c ng + 7c ,fi2) + fi 2 (-2tt 2 + tt 2 + f7 2 )f 



(19) 



(20) 



From these formulae, in the symmetric case that fit = Q c 
and 7c = jc', »[X (1) (0)] « (2 lclab /nl){\ + 2( 7ab7c + 
2tt 2 )/^ 2 ) while dm (1) ]/dA p » -(2 7a6 /r» 2 )(l - 
16 7a 67c/^ 2 t )i which is identical to standard EIT to low- 
est order. However, slightly off resonance in the vicinity 
of the narrow resonances where the absorption is still neg- 
ligible, the affect of these resonances on the dispersion is 
still very significant. 

In Figs. 7 and 8 we consider the effect of the RF Rabi 
frequencies on the group velocity in the vicinity of the 
narrow features at A p ss by plotting the ratio of the 
group velocity for fi;,,f2 c ^ to the group velocity for 
standard EIT. Figure 7 shows the group velocity near 
A p = for zero decoherence and small RF Rabi fre- 
quencies, Of, = Q c = 0.0001 7a fc = IKHz, which leads to 
an almost 100 fold reduction in the group velocity in a 
window of about 50Hz where the absorption is negligi- 
ble. Figure 8 shows the reduction of the group velocity 
and absorption for larger RF Rabi frequencies and finite 
ground state decoherence. As one can see, even in the 
presence of decoherence, one can achieve a reduction in 
the group velocity close to a factor of 10 in a frequency 
window of ~ 5QKHz for the parameters in Fig. 8, while 
at the same time having relatively small absorption. At 
the same time, the dispersion is highly nonlinear implying 
that any pulse propagating near the narrow resonances 
could experience significant reshaping. 



IV. DISCUSSION AND CONCLUSION 

One important aspect that has been ignored in our 
discussion is that of Doppler broadening. To this end, 
we note that the current analysis is limited to that of 
a Doppler free geometry involving co-propagating probe 
and control lasers [2j. Doppler shifts will be unimpor- 
tant for RF/microwave fields. Another alternative to 
avoid Doppler broadening would be to use laser cooled 
atoms in a optical trap [22J. Additionally, we propose 
here that the following energy levels of 87 Rb could be 
used to implement the 5-level geometry (the same tran- 
sitions but with different electronic number could also 
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be used for 23 Na), \a) 
\5S 1/2 ,F = 2, my = 2), \b') = |5S 1/2 
\c) = \5S 1/2 ,F = 2, my = 0), and \c') = 
1, ray = 0, ±1) where the probe laser has a 
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for very small RF Rabi frequencies Clt 
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As in the previous graph, the group velocity is measured rel- 
ative to the group velocity for fij, = fi c = but with all other 
parameters being the same, which we denote as veit- 



while the control laser is a + polarized. We note that the 
five level system could also be implemented using tran- 
sitions between the F = 1 — > F' = manifolds, namely: 
|o) = \5P 3/2 ,F = 0,ray = 0), |6) = \5S 1/2 ,F = l,m f = 
1), \b') = |55 1/2 ,F = l,ray = 0), \c) = \5S 1/2 ,F = 
1 )TO/ = -1), and \d) = \5S 1/2l F = 2, ray = -2,-1,0) 
again with o~ polarized probe and a a + polarized con- 
trol laser. We note that an additional off-resonant laser 
could be used to shift the energy level of \b') relative to 
\b) and \c) via the AC Stark effect so that an RF field 
would only resonantly induce spin flips between \b) and 
\b') and not between \b') and |c). 

In conclusion we have analyzed the linear susceptibil- 
ity in a five level system and showed that the emergence 
of two dark resonances inside of the EIT transparency 
window offers enhanced control over the dispersion and 
absorption spectrum due to the controllable widths and 
separation of the resonances. We have argued that nar- 
rower transparency windows and slower group velocities 
can be achieved using this scheme as compared to stan- 
dard EIT. 

This work is supported in part by the National Science 
Foundation. 
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FIG. 8: (Color Online) (a) The group velocity, v g , in a window 
near A p = for f2 M = 2y a b and Qf, = S7 C = 0.0l7 a and 
~fc — ~fc — 0.000l7a (Blue long dash line); Q b — Q c — 
0.027 a and fc = 7c' = 0.000l7 a (Red short dash line); and 
fib — fl c — 0.0l7a and yc = 7c = (Gold dotted line). In 
each case, the group velocity is measured relative to the group 
velocity for fit = fl c = but with all other parameters being 
the same, which we denote as veit to mean 'standard EIT'. 
We note that the group velocity can readily be related to the 
delay time of a light pulse by r a = L(l/v g — 1/c) where L is 
the thickness of the sample, (b) The imaginary part of the 
susceptibility for the same parameters as in pt. (a). 
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The full statement of the equations of motion deter- 
mining p a B and PaB'i approximated by Eqs. 151 and IT01 
for the case when Af, = A c = 0, and so 9 — 9 C = 7r/4, 
are given here in their general form. We have neglected 
terms proportional to p aa , p cc , p dc >, p ac , p ac >, and p cc >, 
for reasons given in the body of the paper. We find: 



l —flf- = [ a p IT + ~i z 7afc I PaB — [cosObpBB - sm0 b p B , B ) — (cos9 c p C B - sm0 c p C 'B) 



■ dpcB f A A A c A b Of Of 2 . 2 



n, 



■cos9 c paB - i{lc - ic) cos 9 C sin 9 c pcB 



off ocff 



+-£■ sin 6 c paB - i (7C - 7c) cos C sin 9 c p C B 



(21a) 



(21b) 



(21c) 



10 



and 



■ dp a B' f A 



o: 



at 



«7ab ] PaB' + "y ( sin 9bPB>B> - cos d b p B B>) ^- {cos 9 c p C B> - sin 9 c p c >B>) (22a) 



■ dpcs' /. . a c A b nf nf .. 2 „ . 2 \ 

»-^- = ( A ? ~ A ^ + ^ ~ ^ ~ ^ 2~ " * (7C °° c + 7C " Sm c) J /9CB ' 

— — cos9 c p aB > -i(7C -7c) cos 6» c sin 9 c p c >B> 

■ dpcB> ( K A A c A b Of ft? ff v . ,„ 2/ ,A 

Z_ ^T^ "( Ap ~ A ^ + T~T + ^ 2~ ~ * (7C c + 7C ' co c) ) /9c ' s ' 



+-^- su\e c p a B' - i(yc - lc 



2 2 

' - 7c) cos 9 C sin 9 c pcB> ■ 



(22b) 



(22c 



.9/5 



BB 



9t 2 

.dpB'B' hw , ■ 2 



(sin 2 9b(pBB - pB'B>) + sin29bCos29b(pBB' + Pb>b)) 



(sin 9 b (p B B - pB'B') + sin29 b cos29 b (pBB> + PB'b)) 



at 

^'^W' =( ~° ?ff " ^ (3 + cos4 ^))Psb' + ^ ((1 - cos40 b )/5 B , B + sin40 6 (p B , B , - p BS )) 
'^9^ "^ - t -J L{ - 3 + cos4 ^))/5b'B + ^ ((1 - cos4^)/5 BB , + Bm4$ b (p B . B , - ~ P bb)) 



(23a) 
(23b) 
(23c) 
(23d) 
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